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•1—1 , 

^ . We explore gauge actions for lattice QCD, which are constructed such 

^ ! that the occurrence of small plaquette values is strongly suppressed. By 

choosing strong bare gauge couplings we arrive at values for the physical 
lattice spacings of 0(0.1 fm). Such gauge actions tend to confine the Monte 
Carlo history to a single topological sector. This topological stability facili- 
tates the collection of a large set of configurations in a specific sector, which 
is profitable for numerical studies in the e-regime. The suppression of small 
plaquette values is also expected to be favourable for simulations with dy- 
namical quarks. We use a local Hybrid Monte Carlo algorithm to simulate 
such actions, and we present numerical results for the static potential, the 
physical scale, the topological stability and the kernel condition number of 
the overlap Dirac operator. In addition we discuss the question of reflection 
positivity for a class of such gauge actions. 
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1 Introduction 



The gauge action of QCD can be lattice discretised in many ways. One 
requires the naive continuum hmit for all lattice gauge actions to coincide, 
and one usually assumes that then they fall into the same universality class. 
The formulation of the lattice gauge actions is considered less problematic 
than the fermionic part of QCD, in particular due to the absence of additive 
mass renormalisation. Even for the simplest plaquette action, known as the 
Wilson gauge action [1], ^ 

Sw[U] =(3Y,Sp{Up) , Sp{Up) = 1 - ^Re Trf/p , (1.1) 
p 

the lattice artifacts in the scaling behaviour only appear in (9(a^), where 
a is the lattice spacing. This is in contrast to the Wilson fermion, which 
suffers from 0{a) artifacts, as well as additive mass renormalisation, so that 
criticality can only be approached with a tedious fine tuning. 

Nevertheless there are a number of suggestions for improved lattice gauge 
actions. In general the goal is to further suppress the scaling artifacts by 
including more extended closed loops in the discrete formulation of the field 
strength tensor. This is achieved for instance by the tree level improved 
Symanzik gauge action [3], the on-shell improved Liischer-Weisz action [4], 
and by several renormalisation group improved actions [5-7]. 

In the current work we also consider non-standard gauge actions for lattice 
QCD, but we focus on another property, which we want to improve. Our 
intention is to suppress as far as possible the occurrence of small plaquette 
values. Of course, this suppression prevents the gauge field from fluctuating 
as it would for the Wilson gauge action, hence we have to use much stronger 
bare gauge couplings to arrive at a comparable lattice spacing. For practical 
simulations a lattice spacing of 0(0.1 fm) is required; this enables the use 
of a sensible physical volume with lattices of a tractable size. For quenched 
simulations with the Wilson gauge action such a lattice spacing is obtained 
for = G/^fQ ~ 6. We will see that a lattice spacing of the same magnitude 
can still be obtained for the actions that drastically suppress small plaquette 
values, if we drive (3 down to values around and even below 1. 

This observation means that these gauge actions can in fact be used in 
simulations. Now we would like to explain what virtues we expect from it. 

^Up denotes the plaquette variables in a lattice gauge configurations given by the link 
variables Ux,^i G SU{2>). The sum in eq. runs over all plaquettes, see e.g. Ref. [2]. 
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2 Motivations for the suppression of small 
plaquette values 

Small plaquette values Up are expected to be linked to small eigenvalues of 
the Dirac operator, which is most relevant in view of dynamical simulations. 
On the practical side, the suppression of small values for Up should speed 
up dynamical Hybrid Monte Carlo (HMC) simulations, avoid problems with 
the numerical integrator and ease the access to light pions. This would mean 
a further improvement in this direction, in addition to recent algorithmic 
developments [8]. However, these properties are not tested in the current 
work; we leave them for further investigation. 

We do investigate, however, another virtue, which we want to explain 
now. There is a notorious lack of analytical tools to handle QCD at low 
energy, so one often switches to the use of an effective Lagrangian instead. 
In particular, chiral perturbation theory (xPT) is an effective theory in terms 
of the Nambu-Goldstone boson fields of chiral symmetry breaking [9] . This 
approach is very powerful, and it still works if one includes the small explicit 
chiral symmetry breaking due to the finite masses of the light quarks. Then 
one deals with light pseudo Nambu-Goldstone bosons, which are identified 
with the hght mesons. However, by its nature xPT cannot capture all the 
properties of QCD as the fundamental theory. To supplement the missing 
information, in particular the values of the Low Energy Constants in the 
chiral Lagrangian, one has to relate xPT to QCD results, which one tries to 
obtain from lattice simulations, so we need simulations with light mesons. 

In the recent years, a lattice fermion formulation became available, which 
overcomes the conceptual problems related to the chiral symmetry. This 
is the case if a lattice Dirac operator D obeys the Ginsparg- Wilson (GW) 
relation [10, 11]. In its simplest form (and in lattice units, a — 1) the GW 
relation reads 

D^, + ^,D = -D^,D , /.>!. (2.1) 
A* 

The corresponding GW fermions have a lattice modified, but exact chiral 
symmetry at arbitrary lattice spacing, with the full number of generators 
[12]. This property rules out additive mass renormalisation (along with 0{a) 
scaling artifacts). Hence a small (bare) quark mass rriq docs indeed imply a 
small pion mass, oc m^. 
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In practice, there are still difficulties left when we want to approach the 
physical pion mass. To illuminate this point, wc look at the Neuberger 
overlap operator [13], which is an often applied, explicit solution to the GW 
relation. At zero fermion mass it takes the form 



where is the Wilson Dirac operator. Its property -D^ = 75-Dw75 implies 
that Q is Hermitian. A quark mass is added to the overlap operator as 
follows. 



Simulations with this operator are computationally expensive, so that the 
corresponding QCD studies are essentially restricted to the quenched ap- 
proximation up to now. The trouble-maker is the inverse square root, which 
has to be approximated by polynomials of a high degree. The effort to han- 
dle Dfyy is roughly proportional to this degree, which grows — for a fixed 
accuracy of D^v — like the square root of the condition number of Q^. (One 
usually projects out the lowest few modes and treats them separately to lower 
this condition number, but this takes time again.) The lowest eigenvalues of 

are again expected to be linked to the occurrence of low plaquette values. 
We are going to demonstrate an improvement also in this respect for the 
gauge actions that we are going to consider. 

Another obvious problem at small pion mass are the finite size effects. 
In a box of side length L they depend on the product m^L, which should 
be large to keep the finite size effects small, i.e. to be close to the physically 
realistic situation. Then we are in the p-regime, where the p-expansion of 
XPT [14] applies. For instance, for a = 0.1 fm a pion mass of 250 MeV 
corresponds to a correlation length ^ = l/m.,^ of about 8 lattice spacings, 
and to dwarf the finite size effects L should be much larger. Lattices of such 
an extent are very expensive with GW fermions, even quenched, still without 
reaching physical pion masses. 

As a way out one may consider the opposite situation. 



where we are confronted with strong finite size effects. This setting is called 
the e-regime; it allows for the use of a variant of chiral perturbation theory 




(2.2) 




(2.3) 



(2.4) 
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known as the e-expansion [15]. It describes the finite size effects analyti- 
cally, and its formulae can be related to the numerically measured finite size 
effects. The interesting point is that these formulae involve the Low En- 
ergy Constants as they appear in infinite volume. Therefore we can extract 
physically relevant information even from the unphysical e-regime. 

As a peculiarity of the e-regime, the topology plays an important role. 
Observables tend to depend significantly on the topological sector, and pre- 
dictions exist for expectation values in specific sectors [16]. Therefore, for 
numerical measurements it can be advantageous, if not necessary, to disen- 
tangle the topologies, in order to extract maximal information. 

In general, it is not obvious how to define topological sectors on the lattice, 
since all gauge configurations can be continuously deformed into one another. 
A neat definition exists, however, for Ginsparg-Wilson fermions, since they 
have exact zero modes with definite chiralities [11]. Once a GW operator is 
chosen, it has for each gauge configuration a well-defined fermion index 

= — n„ , (2.5) 

where n± is the number of zero modes with positive/negative chirality. In the 
spirit of the Atiyah-Singer Theorem one then uses this index as a definition 
of the topological charge. For independent configurations the distribution of 
V is Gaussian, and its width determines the topological susceptibility. 

Let us assume that we monitor the deformation of a gauge configuration 
to a different topological charge v — defined by the overlap operator given 
in eq. ()2.2j) . In the (real) spectrum of the operator Q such a change of v 
means that some eigenvalue crosses zero (on the topological boundary the 
overlap operator is not defined). For a transition the Monte Carlo history 
has to pass through a region of rather low probability, and in simulations 
one assumes this to happen quickly at some instances in a long history. For 
simulations in the p-regime one tries to sample all topologies and frequent 
transitions are therefore desired. However, in e-regime simulations one would 
often like to collect statistics at one specific value of \v\ to measure an ex- 
pectation value in this sector. For the parameters that have been used in 
the e-regime simulations [17-25], it would be of particular interest to collect 
large sets of configurations with = 1 and 2. The topologically neutral 
sector is problematic due to the frequent appearance of very small Dirac 
eigenvalues, which leads to strong spikes in the Monte Carlo histories of cor- 
relation functions [20] (at |z/| > the non-zero eigenvalues are pushed to 
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higher energies [18,26]). However, a procedure called Low Mode Averaging 
was designed to render also the neutral sector tractable [22] , hence also a cu- 
mulation of configurations with i/ = may be of interest (see also Ref. [24]). 
Charged sectors are required, however, for the method of extracting Low 
Energy Constants solely from the zero mode correlation functions [21,23]. 

A box with K ~ 10 fm^ is suitable, but the width of the Gaussian charge 
distribution is then around iy"^) ^ 10 [19,23], so that most charges vary 
between about —10 and 10. The index measurement by itself is computa- 
tionally expensive, hence identifying a set of, say, (9(1000) configurations in 
one sector is a tedious task — if one uses the standard Wilson gauge action. 

Therefore it is motivated to modify the lattice gauge action such that 
topological transitions are suppressed. ^ The vicinity of a transition (under 
a continuous deformation) is characterised again by the occurrence of small 
plaquette values, hence the technical aim is also in that respect to avoid these 
configurations. This connection was made rigorous first in Ref. [28]. If we 
consider the overlap operator ()2.2p . the square root cannot vanish — and 
therefore topological transitions are excluded under continuous deformations 
— if all the plaquette variables Up in the configurations involved obey the 
inequality (at n = 1) 

where Sp represents the standard gauge action of one plaquette, as specified 
in eq. (jl.l|) . Then the topological structure is continuum-like. 

Later on H. Neuberger showed that this constraint can be relaxed a little 
by increasing the threshold in inequality ()2.6|) to [29] 

1 1 . ^ 

(2.7) 



(1 + 1/V2)rf(rf- 1) 20.5 ■ 

This constraint ensures that the spectrum of (the argument of the square 
root in the overlap formula ()2.2|) ) is strictly positive. ^ In absence of zero 

^Also the use of multi-plaquette gauge actions, such as the those suggested in Refs. [3-7], 
has some impact on the frequency of topological transitions, see e.g. Refs. [27]. But such 
actions are less convenient for a systematic suppression of small plaquette values, so we 
do not consider them here. 

•^The corresponding admissibility condition has also been studied on a non-commutative 
torus in Ref. [30] . 
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modes this condition also guarantees the locahty of the overlap operator (in 
the sense of an exponential decay) [28]. 

On the other hand, the impact of imposing such a restrictive constraint 
strictly could be a severe practical problem. The fluctuations of the gauge 
field would be limited so much that one could only obtain a tiny physical 
lattice spacing, and therefore a tiny physical volume. However, even for 
simulations in the e-regime we have to require that the spatial box length L 
exceeds some lower limit in the range of L > 1.1 fm ... 1.5 fm (depending on 
the exact criterion) [18-24]. 

Here we present numerical experiments with gauge actions which do sup- 
press small plaquette values, but only to the extent that still allows for a 
reasonable physical lattice spacing to be obtained. Then there is no rigorous 
guarantee for topology conservation in the Monte Carlo history. ^ The hope 
is that the transitions are still strongly suppressed, so that the history has 
periods of constant charge, which are sufficiently long to allow us to collect 
many configurations in a specific sector. Moreover, if we can be confident 
that topological transitions rarely happen, most of the index computation 
can be omitted; one would just check after a number of configurations if the 
index has not changed. 

Of course, these configurations should sample independently the observ- 
ables to be measured in a fixed topology. Since we aim at long sequences 
of fixed topological charge, this can also be interpreted as a long topological 
autocorrelation time. Of course, at the same time, we aim at a much shorter 
autocorrelation time for other observables. We repeat that a long topological 
autocorrelation is something one would not wish in general — for instance 
simulations in the p-regime — but in the e-regime this property turns into 
an advantage, if the observables of interest are still weakly autocorrelated. 

*With respect to the topological charge one could object that a Monte Carlo history 
proceeds in discrete jumps, so even with this constraint the charge conservation is not 
absolutely safe. However, this seems like a minor problem, since the charge would still be 
conserved over very long periods in the history, and in simulations of the HMC type the 
few remaining changes could be further suppressed by reducing the step size dr (at higher 
cost, however). 
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3 The gauge actions 



We now describe a number of non-standard lattice gauge actions, which sup- 
press the undesired small plaquette values. As long as the action for very 
smooth configurations — with Sp>0 for all plaquettes — is not altered, the 
naive continuum limit coincides with the one of the Wilson action (jl.ip . and 
with continuum QCD. The suppression becomes strong when Sp reaches the 
value of some parameter e, which one would theoretically choose according 
to eq. ()2.7|) . For practical purposes we will have to relax e to larger values. 
A simple cutoff for Sp at this value would be conceivable, but such a discon- 
tinuity in the action (which would then suddenly jump to infinity) does not 
appear promising. We can still impose a cutoff but let the plaquette action 
diverge continuously as Sp increases towards e, if we modify Sp of eq. 
to the hyperbolic form 



for n > 0. This formulation, with n = 1, was introduced by M. Liischer and 
used for conceptual studies of chiral gauge theories on the lattice [31]. In 
that case, e was of course set to a theoretically stringent value. 

In simulations, this action was first used in the Schwinger model by H. 
Fukaya and T. Onogi [32]. They set e = 1, i.e. far above the theoretical value 
of about 0.29, but they still observed topological stability over hundreds of 
trajectories. ^ 

A theoretical objection against this lattice action was raised by M. Creutz 
[33] . He observed that it does not provide a positive definite transfer matrix. 
Of course, if we rely on the assumption that we are in the same universality 
class as the Wilson action (as the naive continuum limit suggests) then one 
would not worry about that, since the non-positivity would be a cut-off effect 
and we expect positivity to be restored in the continuum limit. Still, a 
problem with this action cannot be strictly excluded. Therefore we kept this 
point in mind in the numerical study and we will comment on it in Subsection 



'^Note that the factor 1/3 in the term for Sp of eq. is actually l/Nc for general 

SU{Nc) or U{Nc) gauge groups. Moreover, the theoretical bound for e in d = 2 is a factor 
6 larger than in c? = 4, as eqs. (|2.6|) and (|2.7() show (this factor is due to a summation 




(3.1) 



4.3. 




in the term for 1/e). 
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The infinite part in action (j3.1|) means that certain steps that the HMC 
algorithm suggests have to be rejected for sure. Therefore we also have to 
verify with special care that the acceptance rate is sufficiently high. 

This motivated us to consider further variants of gauge actions, which 
also suppress the probability of plaquette actions Sp > e, but which do not 
render a violation of this constraint completely impossible. Examples for 
such actions are the "power actions" and the "exponential actions", 

S!°niUp) = SpiUp) + ^Sp{Upr , (3.2) 
StJiUp) = Sp{Up) ■ exp{SpiUpr/e} {n > 0) . (3.3) 

In our numerical studies we included the actions S^^^ and S'^Y- O^^^ 
preliminary results were reported in Refs. [34] and a comprehensive presen- 
tation will be given in the next Section. Further results along these lines for 
quenched QCD can be found in Refs. [35,36]. 



4 Numerical results 

Actions of the types ()3.H) . ()3.2j) and ()3.3|) depend non-linearly on the link vari- 
ables Ux,fj.- As a consequence, the heat-bath algorithm and over- relaxation 
cannot be applied straightforwardly. Therefore we use the local HMC al- 
gorithm, which was introduced in Ref. [37]. Since these actions are still 
composed of separate contributions by the single plaquettes, the force in the 
local HMC algorithm is a simple modification of the corresponding force for 
the Wilson action, 

r,^,, _ ^S^^'jUp) _ ^w... ^ 1 + ^SpjUp) 

SSIJjUp 
where F^(f/, 



= = F^iUp) . (l + ^SpiUpr-') exp {i^p(f/p)"} (4.2) 

SSp{Up) 



is the force of the Wilson action, and its modification by the second factor is 
of order 0{Sp/e) in both cases. This is illustrated in Figure 
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Figure 1: On the left: The ratio between the hyperbolic plaquette action (for 
n = 1) and standard plaquette action, and the corresponding ratio for the 
exponential actions with n = 1 and n = 8 (the latter is that case we studied). 
On the right: the same ratios for the HMC forces, again as a function of Sp. 

In such simulations, it is also of special importance to check that the 
results do not depend on the starting configuration (once we start in the 
desired topological sector). It would be conceivable that the constraint on 
the plaquette values causes also unwanted obstructions. For all the quantities 
to be considered below, it turned out that this was not the case; an example 
is discussed in Subsection 4.6. 

4.1 Plaquette values 

As a first experiment we considered action (j3.H) and searched for the lines of 
a constant mean plaquette action (Sp) on a 4^ lattice, as P and the action 
parameter e are varied. The result is shown in Figure |21 As we decrease e, 
very small values of /3, i.e. strong bare gauge couplings are needed to keep 
{Sp) constant. 

Next we took a look at the statistical distribution of the plaquette actions, 
and we show corresponding histograms in Figure El By decreasing the values 
of e and (3 we can in fact keep (Sp) approximately constant, while drastically 
suppressing the occurrence of very small plaquette actions. 

However, for a systematic approach to identify a line of a constant physical 
scale, we have to proceed to larger lattices and we determine the physical 
scale based on the static potential, since the plaquette value cannot be used 
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Figure 2: The lines of constant plaquette values in the plane spanned by 1/e 
and (3, on a 4^ lattice. 

for this purpose. This issue will be discussed in the next Subsection. 

4.2 The static potential and the physical scale 

A very well established method for setting a scale in pure gauge theory is 
based on the measurement of the static potential at intermediate distances. 
This potential, and the resulting force, are extracted from the Wilson loop 
correlations at sufficiently large time separations, such that excited states 
can be neglected. 

To suppress the effects of short-ranged fluctuations, we applied APE 
smearing [38] on the spatial links before evaluating the Wilson loops. Then 
we evaluated the static potential 



where VF(r, t) is a rectangular Wilson loop with extension r and i in a spatial 
and in the temporal direction (and V{r) is extracted from an asymptotic 
window in t). Here we applied the variational procedure explained in Refs. 
[39,40]. The force F{r) is then obtained from the discrete gradient of V{r). 




(4.3) 
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Figure 3: Examples for the histograms of plaquette actions Sp{Up) with hy- 
perbolic topology conserving actions S^^^' at various parameters, compared to 
the Wilson action, on a 4^ lattice. 
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The aim is to fix the quantity tq = 0.5 fm [41] by tuning the dimensionless 
term 



To this end, we apphed a local interpolation for the force based on the ansatz 



the Wilson gauge action, which we denote by It corresponds to the 

same lattice spacing for the plaquette action (jl.lj) . where we refer to the 
parametrisation formula in Ref. [39]. 

Let us comment now on the possible errors in this evaluation: 

• The systematic error on the interpolation can be estimated by con- 
sidering a third point and observe the relative deviation between the 
two interpolations; in our case this error turned out to be negligible 
compared to the statistical uncertainties. 

• Previous computations [39] revealed that for a box size L > 3.3 vq 
the finite size effects for the computation of can be safely neglected. 
The same work observed that for the Wilson action at /3 = 5.95 and 
L ^ 2.4 ro the finite size artifacts of the force amount to ~ 3%. In our 
study we deal with L = 16 ~ (2.2 . . . 2.4) vq. Hence we assume finite 
size effects for ro to be of this order as well. 

• The errors quoted on Vo/a in Table Q and |21 are purely statistical (they 
were computed by the jackknife method). The same errors are shown 
in Figures 4 and 5, which we comment on below. The extent of these 
errors is acceptable in this context, since a precise determination of ro 
was not the purpose of this work, so we did not aim at high statistics. 

• A way to check the lattice artifacts is to compare the short distance force 
at finite lattice spacing with the one extrapolated to the continuum 
limit in Ref. [40]. In particular we measured the ratio 



r'^F{r)\r=r(c) = c , ro = r(1.65) . 



(4.4) 




A(r/ro) 



r 



^F(r/ro) - r^F(r/ro)|c 
r2F(r/ro)|c 



(4.6) 
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1/e 


(3 


ro/a 


Pw 


(It 


^plaq 


/top 


_plaq . f 
' J top 


acc. rate 


0. 


6.18 


7.14(3) 


6.18 


0.1 


7(1) 


2.2(13) e-2 


^ 1.5 e-1 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.1 


2.2(1) 


3.0(23) e-3 


~ 6.6 e-3 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.05 


2.0(1) 


2.9(11) e-3 


^ 5.8 e-3 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.01 


2.2(1) 


3.5(8) e-3 


^ 7.7 e-3 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.005 


2.3(2) 


2.8(15) e-3 


^ 6.4 e-3 


> 99 % 


1.18 


1. 


7.2(2) 


6.18(2) 


0.1 


1.2(1) 


2.0(12) e-3 


^ 2.4 e-3 


> 99 % 


1.18 


1. 


7.2(2) 


6.18(2) 


0.02/0.01 


1.3(1) 


1.6(7) e-3 


^ 2.1 e-3 


> 99 % 


1.25 


0.8 


7.0(1) 


6.17(1) 


0.1 


1.1(1) 


2.3(13) e-3 


^ 2.5 e-3 


> 99 % 


1.52 


0.3 


7.3(4) 


6.19(4) 


0.1 


0.8(1) 


9.0(28) e-4 


7.2 e-4 


^ 95 % 


1.64 


0.1 


6.8(3) 


6.15(3) 


0.1 


1.0(1) 


1.3(7) e-3 


^ 1.3 e-3 


^ 65 % 


1.64 


0.1 






0.05 


0.7(1) 


2.3(13) e-3 


^ 1.6 e-3 


^ 78 % 


1.64 


0.1 






0.025 


0.6(1) 


3.5(20) e-3 


^ 2.1 e-3 


^ 93 % 


1.64 


0.1 






0.001 


0.5(1) 


3.7(23) e-3 


^ 1.9 e-3 


^ 99 % 



Table 1: Results for the hyperbolic actions S^^^ , defined in eq. for 
different values of e and (5, on a 16^ lattice. We first show the ratio r^ja, 
which fixes the physical scale. For comparison we also display the (3 -values 
13y^, which leads to the same physical scale for the Wilson action M.l]) [39]. 
The trajectories were all of length 1 and divided into HMC steps of length dr. 
For the plaquette values this leads to a mean autocorrelation time r^^^*^, which 
we show as an example for the autocorrelation of a non-topological quantity. 
The topological stability, on the other hand, is measured by the frequency 
of topological transitions. More precisely, /top is the number of topological 
jumps (determined from cooling), normalised by the number of trajectories. 
Its product with t^^^'^ characterises the dominance of the topological autocor- 
relation. Finally we give the acceptance rate of the local HMC algorithm. 
For each set of parameters in this Table we collected at least 200 thermalised 
configurations spaced by 50 trajectories each. A detailed discussion is given 
in the following Subsections. 
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1/^ 




ro/a 




dr 




/top 


rP'^^ ■ /top 


acc. rate 


500 


0.044 


> 8.5 


> 6.30 


0.015 


0.6(1) 


2.5(36) e-4 


1.5 e-4 


^ 99% 


600 


0.0134 


8.0(2) 


6.26(1) 


0.015 


0.6(1) 


2.5(23) e-4 


1.5 e-4 


^ 99% 


1000 


0.004 


> 9 


> 6.34 


0.03 


0.6(1) 


0(0) 




^ 25% 


1000 


0.00113 


7.9(1) 


6.25(1) 


0.015 


0.6(1) 


1.7(23) e-4 


^1.0 e-4 


« 99% 



Table 2: Results for the exponential actions , defined in eq. \3. /or 
different values of e and /3, on a 16^ lattice. As in Table 1 we first show 
the ratio r^/a and the Wilson (3 -value /5w; which leads to the same physical 
scale for the Wilson action M.l]) (here the finite size effects in the evaluation 
of Ti^/a may he sizable). For different HMC steps dr we then give the mean 
plaquette autocorrelation time r^^^^, the frequency of topological transitions, 
/top, its product with r^^'^^ and the acceptance rate. The number of measure- 
ments forr^/a was at least 200 in each case. Further comments are added 
in Subsections 4-4 4-6- 



where r^F(r/ro)|c denotes the continuum limit. The results for the 
action are shown in Figure 4. At short distances, lattice artifacts 
are below 15% for all the different values of 1/e and (3 that we included. 
Moreover, one observes that discretisation errors grow substantially for 
increasing values of 1/e, as expected. This indicates that choosing even 
larger values of 1/e — while keeping the physical lattice spacing fixed 
— could introduce sizeable cutoff effects. 

For the Wilson action it turned out that the lattice artifacts can be re- 
duced considerably by using a so-called "tree level improved" definition 
of the force. For this purpose, one defines an improved distance rj such 
that the force does not contain lattice artifacts at tree level, 

If we adapt this method, the procedure described before leads to the 
results shown in Figure 5. For the largest values of 1/e one observes 
some reduction of the lattice artifacts, whereas they seem to increase 
for the smallest 1/e. This behaviour is not totally surprising, since it 
has been observed in other cases that this improvement is not always 
guaranteed [42]. 
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Figure 4: Lattice artifacts for the action S^^i . The plot shows the relative 
deviation of the force r'^F{r) from the continuum results, given in eq. (|^.6| ). 



We also observe that the use of tree-level improved observables does 
not change the results for ro/a itself (within the statistical errors). 

• For a comparison of the scaling quality one may, for instance, refer to 
the Iwasaki action [5] (at ro/a ^ 6.0) and the DBW2 action [6] (at 
ro/a ~ 5.5) at a distance r/rg ~ 0.3: in these cases the lattice artifacts 
were found to be of order ~ 10% [42]. Those discretisation errors are 
therefore comparable to our results for the actions S^J^. 

We repeat that our statistics is modest, and our intention in this analysis 
was only to check whether the errors and lattice artifacts are reasonably 
under control. This can be confirmed from the results shown Figures 4 and 
5, and the accuracy is fully sufficient for our purposes. 



4.3 The issue of reflection positivity 

As we anticipated in Section 3, Ref. [33] pointed out that actions of the 
type (|3.1|) do not provide a positive transfer matrix, as it was observed [43] 
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Figure 5: Lattice artifacts for the action S^^^ with a tree level improved 
definition of the force \4- The plot shows the relative deviation of the 
force r'^F{r) from the continuum results. 
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a long time ago for the so-called Manton action [44]. The positivity of the 
transfer matrix is assured if positivity holds for both, "site reflection" and 
"link reflection" (e.g. reflections on the hyperplanes X4 = and = 1/2). 

However, this still leaves the possibility of obtaining a positive squared 
transfer matrix. In fact, this is the case if the action obeys at least one of 
the two reflection positivities [2]. We have checked for the actions involved 
in our study, given in eqs. (j3.1|) and (j3.3|) . that site reflection positivity does 
in fact hold. 

From the practical point of view, the lack of positivity can be reflected 
in an irregular behaviour of the effective potential V{r, t) at short time sepa- 
ration, as observed in Ref. [41] for actions involving rectangular loops (like 
those suggested in Refs. [4-6]). This is a lattice artifact which does, in 
principle, not constitute a problem, as long as one keeps far enough from the 
cutoff. ^ In the case of the actions considered in this work and with our 
statistical precision, such an irregular behaviour could not been observed in 
the extraction of the static potential. 

4.4 Acceptance rate 

Also the problem related to the acceptance rate has been mentioned in Sec- 
tion 3. This point motivated us to consider also a set of exponential actions 
of the type (j3.3|) . in addition to the hyperbolic actions, and the corresponding 
results are given in Table |21 In both cases, the acceptance rate is very high 
for most of the actions we studied. It drops, however, if one pushes for very 
low values of e (along with an extremely small (3). As the last lines in Tabled 
show, the acceptance rate can actually be driven up again even at 1/e = 1.64 
by using very short HMC steps. However, this cannot be considered a so- 
lution, because it increases the costs (especially in the dynamical case), and 
also the frequency of topological transitions rises again (c.f. Subsection 4.6). 

Therefore, this property sets another limit on the suppression of the small 
plaquette values, in addition to the scaling of the static potential at short 
distances. 

^Problems can arise in the application of the variational method, where one has to 
choose a small reference time. 
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1.5 


1 


453(101) 


360(30) 


325(11) 


294(7) 


1 


1.18 


390(34) 


328(18) 


302(10) 


281(9) 


0.8 


1.25 


439(89) 


341(23) 


311(15) 


285(9) 


0.3 


1.52 


369(41) 


301(13) 


280(7) 


263(5) 


0.1 


1.64 


433(82) 


342(18) 


315(14) 


293(8) 



Table 3: Condition numbers Cn of the operator in the square root of 
the Neuberger overlap operator at n = 1.6, after projecting out the leading 
n — 1 modes of . In this comparison we always considered configurations 
generated by the local HMC algorithm with dr = 0.1. 

4.5 The condition number for the overlap operator 

This Subsection discusses the condition number of the operator Q^, which is 
crucial for the computational effort required to deal with the overlap operator 
in eq. ()2.2|) . Table 01 collects our results for the action Si^f compared to S'w, 
at the values of e and (3 corresponding to approximately constant physics 
according to Table ^ Similar results were presented in Refs. [36] , which also 
include first tests with dynamical overlap fermions. 

In our study we varied the parameter /i by units of 0.1 and found the 
optimal condition numbers for all gauge actions involved at /z = 1.6. For the 
case of 10 eigenmodes of projected out, this property can be seen in the 
upper plot of Figure El Hence we compare the condition numbers at fi = 1.6 
for different gauge actions in Table El and in the lower plot of Figure El These 
results are based on 30 configurations in each case. 

We show the condition numbers 

Cn:=Amax/A„, (Amax, A„ : largest rcsp. u*"" eigenvalue of Q^) (4.8) 

which are relevant after projecting out the leading n — 1 modes of Q^. ^ 
We see that the c„ are indeed lowered as 1/e increases, which reduces the 

^We have checked that the polynomial degree for a fixed precision of the overlap oper- 
ator is oc (to a high precision). At the side-line, we also observed that there does not 
seem to be to any dependence of the condition numbers c„ on the topological sector. 
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Condition number c^^ at different values of |j. 
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Figure 6: Comparison of the condition number Cu for various values of n (on 
top), and Cn for running n at fi = 1.6 (below), for different types of gauge 
actions. (The parameters /i and Cn are defined in eqs. and Ii4.l^ .) 
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effort for overlap fermion simulations. * If n is around 20, then — for the 
smooth configurations that we considered here — the gain compared to S'w 
is only moderate. However, for the hyperbolic actions S^^f the number of 
these modes can be reduced drastically without much loss in the condition 
number of the remaining operator. This is in contrast to Sw, and it matters 
in applications, since the special treatment of each of these projected modes 
also takes computation time (although this is typically a minor part of the 
total computational effort). 

4.6 Topological charge stability 

For a quick analysis, we used the cooling method [46] to estimate the topo- 
logical charge. The resulting topological stability over the trajectories is 
included in Tables [T] and |21 In independent tests we evaluated for a subset 
of the configurations the overlap indices setting fi = 1.3, 1.4, 1.5, 1.6 and 
1.7. Since we are dealing with smooth configurations, it does not come as a 
surprise that we found an excellent agreement of more than 98 % for all these 
definitions of the topological charge, i.e. the charges obtained with cooling 
and the overlap index at any of the parameters fi hsted up above. ^ Hence 
the results in the Tables are relevant for the overlap index as well. 

As a more direct illustration, we show typical histories of the topological 
charge (defined by cooling) for different parameters in Figure |7| To measure 
the stability of the topological sector, we monitored the number of charge 
changes normalised by the number of trajectories. We denote this parameter 
as the "frequency of topological transitions", /top, and give results in Tables 
□ and 121 

As we already mentioned in Subsection 4.4, for a fixed dr the acceptance 
rate drops when we increase 1/e up to 1.64. Although this problem can be 
alleviated by reducing dr, it indicates that at this point the system tends to 
run too often into forbidden regions, beyond the admissibility cutoff. There- 
fore we also explored an action, which does not have any strictly forbidden 
plaquette values. For instance, for the actions in eqs. (j3.2p and (j3.3p the 
suppression of small plaquette values rises smoothly. Our results for the ex- 
ponential action S'^'^q are collected in Table |21 We see that these actions 

^Alternatively, lower condition numbers c„ can also be achieved by inserting an im- 
proved kernel Q into the overlap formula, see for instance Refs. [23,45]. 

^On the other hand, if we decrease /3w for instance to 5.85, the charge Qtop depends 
significantly on the method of its determination [23] . 
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Figure 7: Typical histories of the topological charge for the actions S'w ( on 
top) and S^^f at various combinations of e and (3, as in TableU[ We show 
results obtained on a 16^ lattice with the HMC step size dr = 0.1. We see 
that an increased 1/e keeps the charge more and more stable. 
The charge was measured by cooling once in 50 trajectories, except for the 
plot at (3 = 6.18, where measurements were made every single or every 5 
trajectories. Notice, on the other hand, that in the plot at (3 = 6.19 the as- 
sumption that the separation of measurements is much larger than the typical 
distance of topology changes is not justified, and the frequency of topological 
transitions based on this plot would be underestimated. 
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do allow us to render the topology somewhat more stable, and again small 
HMC steps help us to keep the acceptance rate high. However, it is difficult 
to fulfil the requirement ro/a<7, although we already chose extremely low 
values for (3. Therefore we did not push further into that direction. 

Let us add some technical aspects about the evaluation of /top- Although 
measuring the charges by cooling is rather cheap, it could still not be eval- 
uated on each trajectory (since we were dealing with quite long histories). 
A reliable determination of /top can only be done if the number of trajec- 
tories, which are skipped between two measurements of the charge Qtop, is 
much less than the typical number of trajectories over which Qtop remains 
constant. It turned out that for Sy/ it was sufficient to cool one configuration 
in 5 trajectories. For the gauge actions at 1/e > 0, one configuration out of 
50 trajectories was typically enough. 

The error on /top was estimated only in a crude way. This is done by 
counting the transitions in 5 sub-histories and taking the standard deviation 
from these 5 samples. The idea is inspired by the jack-knife method, but the 
difference is that the sub-histories have to consist of contiguous elements. 

Of course, we also need to consider the effect of 1/e > on the autocorre- 
lation time with respect to non-topological quantities. One could be worried 
that an improved topological stability comes along with a longer autocorre- 
lation for other observables as well. Our consideration of the plaquette value 
indicates the opposite: its autocorrelation time decreases significantly with 
increasing 1/e, see Tables [T] and |21 

Another conceivable problem could be bad ergocity properties even within 
one topological sector as l/e is switched on. We checked this by perform- 
ing simulations from independent starting points and found that the mean 
plaquette values agree to a very high precision, see Table El and Figure |H1 

Of course, such tests should be extended to further observables, but our 
results for the plaquette value are encouraging. 

5 Conclusions 

The conservation of the topological charge can be implemented in the lattice 
gauge action to some extent. There are various obstacles preventing a strict 
implementation, such as scaling artifacts and the acceptance rate. Still, we 
succeeded in obtaining long sequences of a stable topological charge, although 
it cannot be fixed strictly, as expected. Our comparison to the behaviour of 
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Table 4: Comparison of mean plaquette values (Up) for different parameters 
and different starting points (out of 10 000 trajectories in a volume 16"^). 
The decreased plaquette autocorrelations lead to a much more precise deter- 
mination of (Up). We see a remarkable agreement up to this high precision 
for different starting points, even in different topological sectors. 
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Figure 8: Comparison of short portions of plaquette histories for different 
combinations of j3 ande. 
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the plaquette value suggests that the topological autocorrelation time exceeds 
by far the autocorrelation of other observables. This property facilitates the 
collection of many configurations in a specific topological sector, for actions 
which are perfectly acceptable. We also observed that such actions do not 
seem to suffer from any conceptual problems. In particular the resulting 
static potential is fully consistent with the right continuum limit, and it 
has relatively mild lattice artifacts. We also showed that at least a positive 
squared transfer matrix can be defined. 

Moreover we found benefits of this action for the kernel condition num- 
ber of the overlap operator, which allows for a somewhat faster evaluation 
with a fixed accuracy. We finally remark that our findings are in accordance 
with recent results of Refs. [36]. Further virtues, in particular in view of the 
simulations with dynamical quarks, still remain to be explored. 
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